ABSTRACT
Background: COVID-19 waves caused by specific SARS-CoV-2 variants have occurred globally at different times. We focused on Omicron variants to understand the genomic diversity and phylogenetic relatedness of SARS-CoV-2 strains in various regions of Pakistan. Methods: We studied 276,525 COVID-19 cases and 1,041 genomes sequenced from December 2021 to August 2022. Sequences were analyzed and visualized using phylogenetic trees. Results: The highest case numbers and deaths were recorded in Sindh and Punjab, the most populous provinces in Pakistan. Omicron variants comprised 95% of all genomes, with BA.2 (34.2%) and BA.5 (44.6%) predominating. The first Omicron wave was associated with the sequential identification of BA.1 in Sindh, then Islamabad Capital Territory, Punjab, Khyber Pakhtunkhwa (KP), Azad Jammu Kashmir (AJK), Gilgit-Baltistan (GB) and Balochistan. Phylogenetic analysis revealed Sindh to be the source of BA.1 and BA.2 introductions into Punjab and Balochistan during early 2022. BA.4 was first introduced in AJK and BA.5 in Punjab. Most recent common ancestor (MRCA) analysis revealed relatedness between the earliest BA.1 genome from Sindh with Balochistan, AJK, Punjab and ICT, and that of first BA.1 from Punjab with strains from KPK and GB. Conclusions: Phylogenetic analysis provides insights into the introduction and transmission dynamics of the Omicron variant in Pakistan, identifying Sindh as a hotspot for viral dissemination. Such data linked with public health efforts can help limit surges of new infections.
Subject(s)
COVID-19ABSTRACT
ObjectivesPakistan has a high infectious disease burden with about 265,000 reported cases of COVID-19. We investigated the genomic diversity of SARS-CoV-2 strains and present the first data on viruses circulating in the country. MethodsWe performed whole-genome sequencing and data analysis of SARS-CoV-2 eleven strains isolated in March and May. ResultsStrains from travelers clustered with those from China, Saudi Arabia, India, USA and Australia. Five of eight SARS-CoV-2 strains were GH clade with Spike glycoprotein D614G, Ns3 gene Q57H, and RNA dependent RNA polymerase (RdRp) P4715L mutations. Two were S (ORF8 L84S and N S202N) and three were L clade and one was an I clade strain. One GH and one L strain each displayed Orf1ab L3606F indicating further evolutionary transitions. ConclusionsThis data reveals SARS-CoV-2 strains of L, G, S and I have been circulating in Pakistan from March, at the start of the pandemic. It indicates viral diversity regarding infection in this populous region. Continuing molecular genomic surveillance of SARS-CoV-2 in the context of disease severity will be important to understand virus transmission patterns and host related determinants of COVID-19 in Pakistan.